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We describe a mean field technique for quantum string (or dimer) models. Unlike traditional mean 
field approaches, the method is general enough to include string condensed phases in addition to the 
usual symmetry breaking phases. Thus, it can be used to study phases and phases transitions beyond 
Landau's symmetry breaking paradigm. We demonstrate the technique with a simple example: the 
spin-1 XXZ model on the Kagome lattice. The mean field calculation predicts a number of phases 
and phase transitions, including a z = 2 deconfined quantum critical point. 
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I. INTRODUCTION 

Recently, several frustrated spin systems have been dis- 
covered with the unusual property that their collective 
excitations are described by Maxwell's equations. [1-6] 
These light-like collective modes can be traced to the 
highly entangled nature of the ground state. In these 
systems, the low energy degrees of freedom are not in- 
dividual spins, but rather string-like loops of spins. The 
ground state is a coherent superposition of many such 
string-like configurations - a "string condensate." It is 
this "string condensation" in the ground state that is 
responsible for the emergent photon - just as particle 
condensation is responsible for the phonon modes in a 
superfluid. [3, 7-9] 

While this qualitative picture is relatively clear, quan- 
titative results on string condensation and artificial light 
are lacking. The above models have only been analyzed 
in limiting and unrealistic cases. We do not know any 
realistic systems with emergent photons. The problem 
is that we are missing a good mean field approach for 
exotic states. Current mean field theory approaches can 
only be applied to symmetry breaking states with local 
order parameters. They are useless for understanding 
string condensed states which are highly entangled and 
have nothing to do with symmetry breaking. 

In this paper, we address this problem. We describe a 
mean field approach that can be applied to both symme- 
try breaking and string condensed states. We hope that 
this technique can be used to identify conditions under 
which string condensation may occur, and to help further 
the experimental search for emergent photons and new 
states of matter. 

In practice, our approach can be thought of as a mean 
field technique for quantum string (or dimer) models. 
This technique can be used to estimate the phase dia- 
gram of string (or dimer) models, to find the low en- 
ergy dynamics of the different phases, and to analyze the 
phase transitions. It can be applied to any quantum spin 
system with the property that its low energy degrees of 
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freedom are strings or dimers. This includes all the frus- 
trated spin systems cited above. 

We demonstrate the technique with a simple example: 
a spin-1 XXZ model on the Kagome lattice [3]: 

/ (u) 

-J X yJ2( S I S J + Sy i S j) C 1 ) 
(IJ) 

Here I and J label the sites of the Kagome lattice, and 
J2(u) s u ms over all nearest neighbor sites. This model 
provides a good testing ground for the method since the 
low energy dynamics of H is described by a string model 
in the regime J 2 3> J xy ~> | Ji — Ji \ ■ 

The mean field calculation predicts a number of in- 
teresting phases including string condensed phases with 
emergent photons. The string condensed phases are ul- 
timately destroyed once instanton fluctuations are in- 
cluded, but several phases and phase transitions remain 
- including a deconfined quantum critical point. 

The mean field phase diagram for (1) is shown in Fig. 

la. Here, J= ^ + ^P+3(Ji-J 2 ), and g = ^p. For 

large positive J/g, the system is in a paramagnetic phase 
with no broken symmetries, while for large negative J/g, 
the system is in a plaquette ordered phase with broken 
lattice and spin symmetries. The critical point is in the 
universality class of the 3D XY model. This is similar to, 
but slightly different from the phase diagram obtained in 
Rcf. [10]. In that paper, the authors also predict a pla- 
quette phase at large negative J/g, but their candidate 
phase is a resonating plaquette phase which has differ- 
ent symmetries from the frozen plaquette phase shown 
above. 

We also study the model (1) with an additional sec- 
ond nearest neighbor interaction J3 X^(ij)) S^Sj, J3/ 1 g 
0.17. We find a different phase diagram (Fig. lb). For 
large positive J/g, the system is in a paramagnetic phase, 
while for large negative J / g the system is in a stripe or- 
dered phase with broken rotational and spin symmetry. 
The mean field calculation predicts that the phase tran- 
sition is a deconfined quantum critical point described 
by a U(l) gauge theory with dynamical exponent z — 2. 
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FIG. 1: The mean field phase diagram for (a) the XXZ model 
(1) and (b) the XXZ model with additional next nearest 
neighbor interactions. The filled circles denote spins with 
(S z ) 0. The sign of (S z ) alternates around each plaquette 
in the plaquette ordered phase and along each stripe in the 
stripe ordered phase. 




FIG. 2: When Ji - J 2 = J xy = 0, the Kagome XXZ model (1) 
has an extensive ground state degeneracy. All states satisfying 
Stot,A — X^jgA Sj — are ground states. If we view the sites 
of the Kagome lattice as links of a honeycomb lattice, then 
the ground states are collections of closed loops of alternating 
S z = ±1 on the honeycomb lattice. 



However, we cannot rule out the possibility of a first or- 
der phase transition - for either of the two models. 

The paper is organized as follows. In section II we 
describe the string picture for the Kagome model (1). 
In section III we present the mean field approach and 
derive the mean field phase diagram for (1). In section 
IV, we derive the low energy dynamics in each of the 
phases, and in section V we analyze the phase transitions 
in the model. The details of the mean field calculation 
are presented in the Appendix. 



II. STRING PICTURE 

A. Effective string model 

We will study the XXZ model (1) in the regime J2 ^> 
Jxy 3> \J\ — J2I. In that case the low energy dynamics 




FIG. 3: The low energy effective string Hamiltonian (3) is 
made up of two terms: an operator riij which is the string 
occupation number on the link ij and an operator B p which 
creates (or moves) strings along the boundary of the plaquette 
P- 

of H is described by a string model - a close cousin of 
a quantum dimer model. [3] To see this, note that the 
Hamiltonian can be rewritten as 

h = y E^tot.A) 2 + (Ji - J2) E(sf) 2 

A J 

(IJ) 

where S* ot A = X)/eA and A runs over the triangles 
in the Kagome lattice (see Fig. 2). 

Suppose that J xy = J\ — J2 = 0. In that case, H has 
an extensive ground state degeneracy: every state satis- 
fying S^ ot A = for all triangles A, is a ground state. 
To describe these ground states, it is useful to view the 
sites I of the Kagome lattice as the links of a honey- 
comb lattice - whose vertices we will label by i. One 
ground state is the state with Sf = for all links I of 
the honeycomb lattice. Another ground state can be ob- 
tained by alternately increasing and decreasing Sf along 
a closed loop on the honeycomb lattice. In general all the 
ground states are of this form: they consist of collections 
of closed loops of alternating Sf = ±1 superimposed on 
a background of Sf = 0. (see Fig. 2). 

These states can be thought of as configurations of ori- 
ented strings on the honeycomb lattice. To do this pre- 
cisely, we pick an A and B sublattice of the honeycomb 
lattice. For any link I — (ij), we say that I contains an 
oriented string pointing from i £ A to j e B if Sf = +1, 
and from j to i if Sf — — 1 . We say the link is empty if 
Sf = (see Fig. 3). Then the ground states described 
above are in exact correspondence with configurations of 
oriented closed strings on the honeycomb lattice. 

Now consider the case where J xy , J\ — J2 are small but 
nonzero. These terms will split the extensive degeneracy 
described above. The splitting can be described in degen- 
erate perturbation theory by a low energy effective string 
Hamiltonian H e ff. Working to third order in J xy / J2 and 
assuming J xy \Ji — J%\ we find: 

ij P 
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where J = '—f^- + -jf - + 3(Ji - J 2 ), and g = 5p. 
Here, the two operators, riij } B p are operators that act 
on oriented string states. The operator ny is the string 
occupation number on the link ij. That is, for any string 
state \X), riij\X) = \X) if there is a string on (ij) with 
orientation i — > j and riij\X) = otherwise. The oper- 
ator B p acts on the six links along the boundary of the 
plaquette p. It's action is given by 



a a+1 




where a = denotes the empty link and a = +1 
(— 1) denotes the link with a string oriented clock- 
wise (counterclockwise) with respect to p (see Fig. 3). 
If any of a,..., f = +1, B p annihilates the state: 
B p \a,b,c,d,e, f) = 0. One can think of B p as being 
analogous to the dimer flip operators in quantum dimcr 
models. 

In the remainder of the paper, we will focus entirely on 
this low energy effective string model (3). While our dis- 
cussion will focus on strings we would like to emphasize 
that the results tell us about the physics of the simple 
spin model (1). 




FIG. 4: One possible phase diagram for the string model (3). 
The string orientations have been omitted for clarity. 

Clearly, even the qualitative phase diagram for the 
string model (3) is not obvious. There are many po- 
tentially competing phases with very different properties. 
We would like to have a method for determining which of 
these phases are actually realized and for what param- 
eters. It would be particularly useful to know in what 
region the string condensed phase occurs, if at all - since 
this phase may contain gapless photon-like modes. The 
mean field technique described below accomplishes this 
task. 



III. MEAN FIELD APPROACH 



B. Possible phases 

The two terms in the Hamiltonian (3) have a sim- 
ple physical interpretation in the string language. The 
first term, J/3 • ny is a string tension which penalizes 
strings for being long. The second term, g/2J2 p (B p + 
h.c) is a string kinetic energy term. It makes the strings 
fluctuate and gives them dynamics. 

We are interested in the phase diagram of this model. 
We will focus on the case g > for simplicity. There are 
three basic regimes to consider: J/g » 1, |J|/g <C 1 and 
-J/g » 1. 

When J/g 3> 1 the string tension term dominates. In 
this regime, the physics of (3) is clear. We expect the 
ground state to be the vacuum state with a few small 
strings; the system is in a "small string" phase. On the 
other hand, when \J\/g <C 1, the string kinetic energy 
dominates. In this regime, the behavior of (3) is less 
straightforward. One plausible scenario is that kinetic 
energy term favors a ground state which is a superposi- 
tion of many large strings - a "string condensed" state, 
(see Fig. 4) However, it is also possible that the ener- 
getics favor a string crystal phase, or even a small string 
phase. In the final regime, when — J/g ^> 1, the negative 
string tension favors "fully packed" string states - states 
where every point on the honeycomb lattice is contained 
in a string. Again, the physics is not clear. It is possible 
that the system enters a fully packed string crystal phase 
- but it could equally well realize a fully packed string 
liquid state. 



In this section, we describe a mean field technique 
for quantitatively computing the phase diagram for the 
string model (3). As we will see later, it can also be used 
to derive a low energy effective Lagrangian for each of the 
phases, and to analyze the phase transitions. We would 
like to mention that similar approaches have been used 
to study lattice gauge theory. [11, 12] 

A. Variational states 

Our approach is variational - we define a large class 
of variational string wave functions, and then minimize 
their energy (H). Let us begin by describing the varia- 
tional wave functions. The wave functions have a large 
number of variational parameters {zij} indexed by the 
oriented links ij of the honeycomb lattice. For each set 
of {zij}, the corresponding wave function ^{ z } is defined 

by 

ij 

where ny is the occupation number of the oriented link 
ij in the oriented string configuration X. We can see 
that Zij,zji are string fugacitics on the link (ij) for the 
two different string orientations. 

The above variational wave function (4) can accommo- 
date many different kinds of states - including both string 
crystals and string liquids. If zy is periodic, ^{ z } is a 



symmetry breaking string crystal state; if Zij is constant 
for all ij, then ^i z \ is a string liquid state. The varia- 
tional states ^{ z } can even access the two types of string 
liquids described in the previous section - small string 
states and string condensed states - and can capture the 
distinction between them. 

To see this, consider a string liquid state with zij — a. 
The properties of the state ^ a can be deduced from 
the properties of the classical loop gas with statistical 
weight P{X) = \^ a {X)\ 2 = a 2L ^ (here L(X) is the 
total length of all the loops in X). This loop gas is 
the classical 0(2) loop gas on the honeycomb lattice 
and has been solved exactly. It is known that the loop 
gas has two phases separated by a phase transition at 
a c = 2~ 1/4 w 0.84. [13] When a < a c , the classical loop 
gas is in a "small loop" phase, where the typical loop size 
is some finite length scale £. On the other hand, when 
a > a c , the loop gas enters a phase with large loops 
of arbitrarily large size. Intuitively, this means that the 
states with a < a c should be regarded as small string 
states, while the states with a > a c should be regarded 
as string condensed states. 

In the same way it is not hard to see that the wave 
functions ^{ z } can even accommodate simultaneous sym- 
metry breaking order and string condensation. Consider, 
for example, the case where zy is large and nearly con- 
stant, but has a small periodic position dependence. In 
that case, ^{ z } will exhibit both translational symmetry 
breaking and string condensation. 



B. Defining string condensed states 

In the preceding section, we relied on an intuitive pic- 
ture of string condensed states. Here we make our lan- 
guage more precise - formally defining which variational 
states ^{ z } we regard as string condensed. 

To state our definition, we place the wave function 
^{ z } on a thermodynamically large cylinder and consider 
the classical loop gas associated with l^^}! 2 . We define 
a quantity p s by 

Ps = (W 2 ) (W) 2 (5) 

where the expectation values are taken with respect to 
the classical loop gas, and W is the winding number, e.g. 
W counts the number of times the loops wind around the 
cylinder. Our definition is that the state ^{ z } is string 
condensed if p s > and not string condensed if p s = 
in the thermodynamic limit. Thus, p s can be thought of 
as a (non-local) order parameter for string condensation. 

The motivation for this definition is twofold. First, 
it captures our intuitive picture of string condensation: 
states with large fluctuating loops (such as ^> a , a > a c ) 
have p s > while other states (such as a < a c ) 
have p s = 0. Second, it agrees with the expected low en- 
ergy physics of string condensed states. We will see that 
states with p s > generically have a gapless linearly dis- 
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FIG. 5: Monte Carlo results for p s as a function of a, obtained 
on 6 x 6, 9 x 9 and 12 x 12 lattices. In the thermodynamic 
limit, p s vanishes for the small string states and then jumps 
to a finite non-zero value for the string condensed states. We 
estimate the critical point a c ~ 0.82 as the value of a where 
the three finite size curves intersect. 



persing photon-like mode (neglecting instanton effects), 
while other states are gapped. 

In addition to clarifying our discussion of string con- 
densation, p s can be used in numerical simulations to 
distinguish string condensed ^{ z } from normal ^{ z }- As 
a demonstration, we have computed p s numerically for 
the liquid states, Zij — a and found that the transition 
from small string states to string condensed states occurs 
at a c w 0.82 (Fig. 5). In this particular case this com- 
putation was unnecessary, since the transition point was 
already known analytically: a c = 2 -1 / 4 w 0.84. How- 
ever, in more complicated cases (such as states with non- 
uniform Zij), such a numerical calculation would be nec- 
essary. 

Finally, we would like to mention a physical interpre- 
tation of p s which will prove useful later. This inter- 
pretation is based on the duality between the classical 
0(2) loop gas and the classical (2D) XY model. [13] 
This duality maps the small loop phase of the classical 
loop gas onto the disordered phase of the XY model, the 
large loop phase onto the (algebraically) ordered phase 
of the XY model, and the transition at a — a c onto the 
Kosterlitz-Thouless transition. Under this duality, the 
quantity p s corresponds to the supcrfluid stiffness in the 
dual XY model. This means that our previous definition 
of string condensed states can be rephrased as: ^{ z } is 
string condensed if and only if the XY model dual to the 
loop gas \^ z \\ is in the ordered phase. 



C. Mean field phase diagram 

The mean field phase diagram can be obtained by min- 
imizing the ground state energy (H) over all states \& { z \, 
and then identifying the quantum phase associated with 
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(b) 



Plaquette order 



Small string 



-0.43 
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String condensed 

Stripe order . Stripe order . Siring condensed . Small siring 



-1.83 



-0.11 



(b) 



Stripe order 



Small string 



-1.83 



FIG. 6: (a) The phase diagram obtained from the mean field 
approach, (b) The phase diagram after instanton effects have 
been included. The string condensed phases are destroyed 
and only two phases remain. 



the minimum energy ^{ z }- Energy expectation values 
can be obtained in a number of ways - we compute them 
numerically using a variational Monte Carlo method. 

We have applied this technique to (3), using the energy 
minimization procedure described in Appendix A. We 
find the mean field phase diagram shown in Fig. 6a. We 
find that when J/g > 0.27, the system is in a small string 
liquid phase. When —0.43 < J/g < 0.27 the system is in 
a string condensed liquid phase. When J/g < —0.43, the 
system enters a phase with simultaneous string conden- 
sation and plaquette order. We have not executed sys- 
tematic numerics beyond this point but we believe that 
when J/g becomes sufficiently large and negative, the 
string condensation is destroyed, and the system enters a 
phase with plaquette order and no string condensation. 
(See Appendix B for details on how these results were 
obtained) . 

As we will show later, the true phase diagram is differ- 
ent from the mean field phase diagram due to instanton 
effects. Once these are taken into account, the two string 
condensed phases are destroyed, and only the small string 
and plaquette ordered phases remain (Fig. 6b). This is 
similar to, but slightly different from the phase diagram 
obtained in Ref. [10]. In that paper, the authors predict 
a resonating plaquette phase (Fig. 8b) at large negative 
J / g, while our mean field phase diagram predicts a frozen 
plaquette phase (Fig. 8a). The two phases have the same 
translational symmetry group, but different symmetries 
under spatial and spin rotations. 

We have also computed the mean field phase dia- 
gram from the model (1) with an additional next nearest 
neighbor interaction J3 Yl/Uj)) 5f &j 1 = 0.17. We 

find that the phase diagram shown in Fig. 7a. When 
J/g > —0.11, the system is in a small string liquid phase, 
while when —1.83 < J/g < —0.11 the system is in a 
string condensed liquid phase. When J/g < —1.83, the 
system enters a phase with simultaneous string conden- 
sation and stripe order (Fig. 9). Instanton effects ulti- 
mately destroy the two string condensed phases leaving 
only the small string phase and the stripe ordered phase 
(Fig. 7b). However, as we will see in section (V), an 
interesting deconfined quantum critical point remains - 
the transition between the small string and stripe ordered 
phases. 



FIG. 7: (a). The mean field phase diagram for the XXZ 
model with a small second nearest neighbor interaction 
JsYltUj)) &i Sj > J3/9 — 0.17. (b). The phase diagram af- 
ter instanton effects have been included. 




(a) (b) 



FIG. 8: The mean field calculation predicts (a) frozen plaque- 
tte order for J/g < —0.43. Previous work [10] has suggested 
(b) resonating plaquette order. The two orders are shown 
above: the thickness of the bonds indicates the size of «y ■ Zji 
while the arrows indicate the relative size of Zij and Zj%. 



IV. LOW ENERGY DYNAMICS 

What are the low energy dynamics in each phase? 
From the general string picture, we expect (oriented) 
string condensed phases to be described by compact (7(1) 
gauge theory. Therefore, we expect that the two string 
condensed phases have a gapless photon-like mode (ne- 
glecting instanton effects), while the other phases are 
gapped. But we would like to understand this more 
quantitatively. We use the mean field method to ac- 
complish this task. We construct a low energy effective 
Lagrangian that describes the dynamics of the string col- 
lective modes. 

The idea is based on the coherent state approach (or 
dynamical variational approach). Suppose that for some 
value of J, g, the variational ground state is ^{z}- It " 1S 
useful to label this state - once it's been properly normal- 




FIG. 9: The addition of a second nearest neighbor interaction 
leads to a stripe ordered phase, depicted above. The thickness 
of the bonds indicates the size of Zij ■ Zji while the arrows 
indicate the relative size of Zij and zji. 
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ized - as While !{;%}} is not the exact ground 

state, our approach is based on the assumption that a 
good approximation to the ground state can be obtained 
by taking linear combinations of states \{z%j + Szij}} for 
small Szij. We also assume that low energy excitations 
can be represented by such states. With this assumption, 
we can use \{z~ij + Szij}) as coherent states and use the 
coherent state path integral to calculate the low energy 
dynamics. In general, the Lagrangian for the coherent 
state path integral is given by 

L(Sz, Sz) = ({z + 6z(t)}\{i± - H)\{z + Sz(t)}) (6) 

where the first piece (z + 5z\i4z\z + 5z) is the Berry phase 
term, and the second piece (z + Sz\H\z + Sz) is the usual 
Hamiltonian evolution term. 

We will focus on the simplest case, where the ground 
state is a string liquid: z^ = a for all ij. It is conve- 
nient to parameterize the fluctuations by Zij + Sz^j — 
ae Eij+tA i:j wnere Eij,Aij are arbitrary real numbers. 
Substituting this into (6) and expanding to quadratic or- 
der in E, A gives a Lagrangian 

L= *^^(bij t kiEijAki — Cij,kiEijEki—dij t kiAijAki) (7) 

ij.kl 

Here, the constants bijM, Cij, k i, dij,ki are given by the 
correlation functions 

hjM = -2((nijn k i) - (nij)(n kl )) (8) 
Cij,ki = 2((nijHn k i) + (HnijUki) -2(H)(nijn kl )) 
dij.ki = 2({mjHn k i) - (HriijUki)) 

evaluated for the string liquid state \I> a . 



A. Gauge structure 

One might expect that by analogy with the harmonic 
oscillator action, one could quantize (7) using states of 
the form |{^4jj}) as a complete orthonormal basis. How- 
ever, this is not quite correct. The problem is with our 
starting point: the coherent states |{2y}) are not all dis- 
tinct. The loop structure of the string states means that 
K^ij'}) — K^ij}) ^ i z } an d differ by a transformation 
of the form z%j = e l ( Xi ~ x ^Zij. 

In terms of our parameterization z^ — ae By+lj4i) ' , this 
means that |{j4jj}) = |{j4jj}) if Aij = Aij + \i — \j. 
Thus, the are a many-to-one labeling of states. 

To obtain a true orthonormal basis, we need to treat each 
equivalence class of {^4ij} as a single state. 

This many-to-one labeling is identical to the U(l) 
gauge structure in U(l) lattice gauge theory. Thus the 
above action (7) should be regarded as a U(l) gauge the- 
ory. 



B. Band structure 

To solve (7), we need to find the normal modes. The 
unit cell contains three links with two possible orien- 
tations, so there are six bands. Three of these bands 
correspond to symmetric configurations, Eij — Eij, 
Aij = Aij and three correspond to antisymmetric con- 
figurations Eij = —Eij, Aij = —Aij. The symmetric 
bands are generically gapped so we will focus primarily 
on the antisymmetric bands. 

Two of the antisymmetric bands are exactly flat with 
Wfe = for all k. These bands correspond to pure gauge 
fluctuations - modes of the form Aij = Aj — Xj . They are 
exactly flat because pure gauge fluctuations don't change 
the physical state at all. These modes don't correspond 
to physical degrees of freedom. 

The only physical low energy degree of freedom is given 
by the third antisymmetric band. This band corresponds 
to transverse modes of A, E. In the limit fc — ^ 0, it is 
given by 

A k , rs = -J=( n .(s-r))A k e ik - r (9) 

E k , rs = -^(„.(,- P ))£fc c ** 

Here n _L k is a unit vector, a is the lattice spacing and 
N is the number of plaquettes in the lattice. 

Substituting these expressions into the Lagrangian (7), 
we find 

L = J2(bkE* k A k -c k \E k \ 2 -d k \A k \ 2 ) (10) 
k 

where 

bk = ^ E b r8 , r , a ,(n-(s-r))(n-(s'-r'))e ik < r '-^ 

rs,r's' 

and Cfc, d k are defined similarly. The dispersion can now 
be easily obtained: u> k = 

C. Gapless photon mode 

The above Lagrangian (10) is a 17(1) lattice gauge the- 
ory. Therefore, one might expect that the system always 
contains a gapless photon mode. 

This is not the case. We will now show that the pho- 
ton mode is gapped in certain phases. Specifically, we 
will show that the mode is gapped in the small string 
phase and gapless in the string condensed phase. This is 
exactly what we expect physically based on the general 
string condensation picture. What is perhaps surprising 
is that this confincmcnt/dcconfincmcnt physics can be 
captured without including the compactness of the U(l) 
gauge field. 

To analyze the low energy excitations, we need to con- 
sider the limit k — > 0. In that limit, d k ~ dk 2 a 2 for some 
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constant d. One way to see this is to use the original defi- 
nition of the Lagrangian L. According to that definition, 
dk is proportional to 

d k oc (ae iAh ^\H\ae lA ^') - (a\H\a) 

But it's not hard to see that 

(ae lAk -™\H\ae iAk -™) - (a\H\a) oc ^(cos(F p ) - 1) 



where F p = A\i + A 23 + A 3i + A i5 + A 56 + A e i is the flux 
through the plaquctte p. Since for small fc, cos(_F p ) — 1 <~ 
k 2 a 2 we conclude that dk ~ dk 2 a 2 for some constant d. 

Because dk — > as fc — > 0, there are potentially gapless 
excitations at fc = 0. The presence or absence of a gap 
depends on the behavior of bk and Ck at small fc. To 
understand this behavior, we make use of the duality 
between the classical loop gas and the XY model. [13] 
Under this duality the quantity n rs — n sr corresponds to 
the boson current (s — r) ■ j(r) (j(r) = V#(r)). This 
means that bk can be identified with the current-current 
correlator, b h ~ ((n-j k )(n -j_ k )). 

There are two regimes to consider. In the string con- 
densed phase, a > a c , the XY model is in the ordered 
phase. In this case, the current-current correlator is of 
the form 



Ofei-fe) 



Wk v - k 2 S^ 



(11) 



This means that bk approaches a nonzero constant value 
b <~ p s as fc — > 0. 

On the other hand, in the small string phase, a < a c , 
the XY model is in the disordered phase. In this case, 
the current-current correlator is of the form 



k^k v - k 2 8^ 

k 2 + i/e 



(12) 



where £ is the correlation length. This means that bk ~ 
k 2 £ 2 as fc — > 0. In a similar way, we can derive the 
behavior of Ck as fc — > 0. Note that 

Crs.tu - d rSitu = 4:({H5 rs 5 tu ) - (H)(5 rs S tu )) (13) 

The Fourier transform of the left hand side is Ck—dk- The 
Fourier transform of the right hand side can be expanded 
in terms of current operators as 



-i-qjqj-kjk} 



nK S'lJk • - ( 14 ) 

<3 



By Wick's theorem, this expression is of the same order 
as b\. It follows that as fc — > 0, Ck ~ dk ~ <ifc 2 a 2 in the 
small string phase, and Cfc approaches a nonzero constant 
value c ~ p 2 in the string condensed phase. 

Putting our expressions for 6fc,Cfc,<ifc together, we can 
compute the dispersion relation ojk as fc — > 0. We find 
that in the small string phase, the fc — * excitations are 



string condensed phase, the excitations are gapless: LUk ~ 
a k* |fc|. This gapless mode is the artificial photon that 
we expected in the string condensed state. The photon 
propagates with a "speed of light" v < 



'cd 



We can apply the same arguments to the other phases 
in the mean field phase diagram. A similar calcula- 
tion predicts gapped excitations in the plaquette ordered 
phase and a gapless photon mode in the phase with si- 
multaneous plaquette order and string condensation. 

Clearly, the presence or absence of a gap is directly re- 
lated to the behavior of the string-string correlations at 
small momenta or large distances. Only in the string con- 
densed phase, where the string-string correlations decay 
algebraically, is a gapless photon mode present. Indeed, 
this connection can be made rigorous using a single mode 
approximation argument similar to that of Ref. [14]. One 
can prove a "Goldstone theorem" for (oriented) string 
condensation which asserts that any string state with al- 
gebraic string-string correlations (S rs St u ) ~ 1/ jr* — t\" 
with v < 2 + d (where d is the spatial dimension) has a 
gapless photon mode. 



D. Instanton effect 

To fully understand the low energy physics of the string 
condensed phase, we need to go back to a real space de- 
scription. We restrict our attention to low energy, long 
wavelength fluctuations of the form (9). For slowly vary- 
ing modes like these, the Lagrangian (10) can be written 
in real space as 

L = bJ2 EijAtj -cJ2E^-dJ2 F 2 P (15) 



13 



^3 



where F p = A\i 



gapped: uj k 



On the other hand, in the 



A 23 + A 34 + A 45 + A 56 + A 61 is the 
flux though the plaquette p. Notice this is precisely the 
Lagrangian for U(l) lattice gauge theory. It gives rise 
to light-like collective excitations with a speed of light 

a\fcd 
U b ■ 

This is exactly what we claimed earlier. However, in 
the preceding discussion we neglected an important ef- 
fect. The above lattice gauge theory is actually a com- 
pact U(l) gauge theory: Aij and Ay + 2ir represent the 
same state. Therefore, the magnetic energy term should 
be changed from F 2 to —cos(F p ). This has dramatic 
consequences for the low energy physics. Due to the non- 
perturbative instanton effect, compact U(l) gauge theory 
is always confining in (2 + 1) dimensions. The photon 

mode obtains a finite gap of order A 

V a e o 

where e\ <~ ^Jc/db 2 is the dimensionless gauge coupling 
and K is a dimensionless constant of order 1. [15] 

Physically, this means that fluctuations, in particular 
instanton fluctuations, modify the mean field phase di- 
agram derived earlier. The instanton fluctuations pre- 
vent the strings from obtaining an infinite correlation 
length, and thus destabilize the string condensed phase. 
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FIG. 10: Three types of phase transitions between string con- 
densed phases and other phases are allowed within mean field 
theory. Two are shown above. The first (a) occurs when the 
photon mode at k = becomes unstable. The second (b) 
occurs when a photon mode at some k ^ becomes unstable. 

By similar reasoning, we expect the fluctuations to de- 
stroy the phase with simultaneous string condensation 
and resonating plaquette order. Thus, once we take in- 
stanton fluctuations into account, all that remains arc 
two phases: the small string phase and the plaquette 
phase (Fig. 6(b)). 

While the instanton fluctuations dramatically alter our 
phase diagram in (2 + 1) dimensions, we expect that a 
similar mean field analysis in (3 + 1) dimensions to be 
much more stable. In (3 + 1) dimensions, the only effect 
of quantum fluctuations (such as monopole fluctuations) 
is to modify the position of the string condensation phase 
transition. Therefore, our mean field analysis may be 
most useful in this context. The example described in 
this paper may be viewed as a warm up for such a three 
dimensional calculation. 



V. NATURE OF THE PHASE TRANSITIONS 

It is natural to wonder what the mean field approach 
can tell us about the phase transitions. In this section, we 
discuss this issue, focusing on phase transitions between 
the string condensed phase and adjoining phases. 

Within mean field theory, there are three basic types 
of transitions that can occur between a string con- 
densed phase and a neighboring phase: (1) confine- 
ment/deconfinement transitions, (2) transitions resulting 
from an instability of a k + 1 photon mode (Fig. 10b), 
and (3) transitions resulting from an instability at k = 
(Fig. 10a). Examples of all three of these transitions 
can be found in the simple spin-1 XXZ model. The third 
type of transition is particularly interesting since it is 
described by a z = 2 deconfincd quantum critical point. 



A. Confinement / deconflnement transition 

An example of a confinement /deconfinemcnt transition 
is given by the critical point separating the small string 
and string condensed phases. According to mean field 
theory, this transition is described by a singularity in 
the variational ground state When J/g is large and 



positive, the energetics favor a small string state with 
a < a c . On the other hand, when J/g is smaller, the 
energetics favor a string condensed state with a > a c . 
The critical point occurs when the energetically favored 
a tunes through the critical value a c . 

The behavior of the low energy excitations near the 
critical point can be derived from the effective La- 
grangian (10). First, consider approaching the critical 
point from the small string side. On the small string 
side, the dispersion is linear in k for k ^> a/^ 2 , and levels 
off to a gap proportional to a 2 /£ 2 when k <C a/£ 2 . As we 
approach the critical point, the linear dispersion extends 
to longer and longer length scales and the gap goes to 
zero - resulting in a gapless photon mode. 

The behavior on the string condensed side of the criti- 
cal point is even simpler. In that case, there is no length 
scale other than the lattice spacing a. For all k <C 1/a 
the dispersion is linear in k. As we approach the transi- 
tion, both b (~ Ps) and c (~ p 2 ) approach finite non-zero 
values. Thus, the transition occurs at a finite photon 
velocity v and a gauge coupling e 2 , <~ 1. 

This mean field picture may capture some of the quali- 
tative features of the small string/string condensed phase 
transitions, but it is almost certainly incorrect when it 
comes to a quantitative description. A major source of 
suspicion is that the mean field phase transition origi- 
nates from a singularity in the variational wave function 
^ a itself, not from the fluctuations about this state. As a 
result the mean field exponents for the (2+1) dimensional 
system come from a critical points in (2 + 0) dimensions. 
It seems unlikely that this is the correct quantitative de- 
scription of the critical point. 

This issue is not specific to two dimensions. If the 
mean field technique is applied to a small string/string 
condensed phase transition in (3 + 1) one finds that the 
mean field exponents come from a critical point in (3 + 
0) dimensions (the 3D XY model). Again, this seems 
incorrect. 

The problem is that when we derived the Lagrangian 
(15) we neglected higher order fluctuations, in particu- 
lar the instanton (or monopole) fluctuations of the U(l) 
gauge field. These fluctuations can completely change 
the mean field phase transition. To include them we have 
to treat the U(l) gauge field as a compact gauge field, 
replacing the F 2 term in (15) by — cos(F p ). 

Taking these fluctuations into account gives a more ac- 
curate picture of the small string/string condensed crit- 
ical point. In (2 + 1) dimensions, instanton fluctuations 
destroy the string condensed phase altogether, and there 
is no phase transition at all. On the other hand, in 
(3 + 1) dimensions, monopole fluctuations give a com- 
pletely different mechanism for a small string/string con- 
densed phase transition - namely monopole condensation. 
This monopole mediated transition (which is known to 
be weakly first order [16, 17]) will likely preempt the (un- 
physical) mean field critical point. 
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FIG. 11: The dispersion u>k of the photon-like mode for dif- 
ferent values of J/g. As J/g approaches the critical point 
(J/q)c — —0.43, the energy of the k = Q photon goes to 0. 
Thus, this critical point is an example of the class of phase 
transitions depicted in Fig. 10b. 



k=0 



FIG. 12: The Brillouin zone for the honeycomb lattice. The 
instability to plaquette order occurs at modes (Bl) with wave 
vectors k — ±Q. 



where £(r),#(r) vary slowly on the scale of the lattice 
spacing. Substituting these expressions into the La- 
grangian (7), and integrating out the 9 field gives an 
action of the form 



(\d t (\ 2 - v 2 \d x (\ 2 ) - A\(f 



(18) 



where A oc cq. 

This is still not quite the correct critical theory. To 
get the full theory, we need to take fluctuations in ( 
into account by including higher order terms. Following 
the analysis preceding (B2) , we include the most general 
terms consistent with the lattice symmetries: 



C 



, 2 \d x (\ 2 )-A\(\ 2 -B\C\ 4 -C(( e +h.c.) 

(19) 

The result is the critical theory for a Zq symmetry break- 
ing transition - not surprising, given the symmetries of 
the two phases. At the transition, the £ 6 term is irrele- 
vant and the critical point is in the universality class of 
the 3D XY model. 

We would like to mention that there is another possibil- 
ity that we cannot rule out: the phase transition could be 
first order. Because we implemented the restricted min- 
imization procedure described in Appendix A we cannot 
resolve this question. However, in principle the mean 
field technique can address this issue. One simply needs 
to use a more general minimization procedure. 



B. Transition via instability at k ^ 

The S = 1 XXZ model also contains an example of a 
transition resulting from an instability of a k ^ pho- 
ton mode (Fig. 10b). This example occurs at the critical 
point separating the string condensed state from the state 
with simultaneous plaquette order and string condensa- 
tion. In the mean field picture, this transition occurs 
when the mode at k = ±Q becomes unstable (Fig. 11, 
12). The critical point is described by tuning one of the 
coefficients cq in (10) through (Fig. 15). When cq > 
the ground state is a string condensed liquid state. When 
cq < 0, nonzero E±q is energetically favorable, and the 
ground state acquires plaquette order (in addition to the 
string condensation). 

Again, we need to include fluctuations to obtain the 
correct physics. Instanton fluctuations gap the photon 
modes, destroying the string condensation on both sides 
of the transition. The result is that the transition is 
actually a simple symmetry breaking transition -between 
a small string state and a plaquette crystal state (Fig. 
6b). 

To obtain the critical theory, note that due to the in- 
stanton fluctuations, the only low lying modes are those 
with k rts ±Q. These modes can be parameterized as 

E rs = ((r)E 

-\-Q,rs 

+ C(r)E_ Q , ra (16) 
A rs = 9(r)A +Qtra + e*(r)A_Q !ra (17) 



C. Transition via instability at k = 

The third type of transition - where a photon mode 
at k = becomes unstable (Fig. 10a) - does not occur 
in the nearest neighbor XXZ spin model (1). However, 
with an additional second nearest neighbor interaction 
J3 J2((ij)) sucn a transition does occur. 

With the addition of the second nearest neighbor term, 
the mean field phase diagram changes to the one shown 
in Fig. 7(a). The phase adjacent to the string condensed 
phase is not plaquette ordered, but instead is stripe or- 
dered (Fig. 9) with simultaneous string condensation. 
The transition occurs when the k — mode 

E rs = n-{s- r) (20) 



becomes unstable (Fig. 13). The mean field critical the- 
ory is described by tuning the coefficient c = lim^^o c fe 
through in the effective Lagrangian (10). 

It is useful to rewrite this theory in real space. Since we 
arc only interested in the low lying modes with k small, 
we can make the approximation bf. w b, Ck ~ c + ek 2 , 
dfc w dk 2 . Going to a real space, continuum description 
gives a Lagrangian of the form 

C = b'EA - c'E 2 - d'{V x A) 2 - e'(V x E) 2 (21) 
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FIG. 13: The dispersion ui h of the photon-like mode for dif- 
ferent values of J/g and Jz/g = 0.17. As J/g approaches the 
critical point (J/g) c = —1.83, the velocity of the photon at 
k — goes to 0. Thus, this critical point is an example of the 
class of phase transitions depicted in Fig. 10a. 



Thus, the z = 2 critical point (21) is potentially sta- 
ble (just as in the previous section, we cannot rule out 
the possibility of a first order phase transition). This 
is an example of a deconfincd quantum critical point. 
While the two adjoining phases differ by simple Z 6 sym- 
metry breaking, the phase transition is not captured by 
a Landau-Ginzburg- Wilson action. Instead, the critical 
theory is described by gauge fluctuations which become 
deconfined only at the critical point. 

This deconfinement is physically reasonable - the crit- 
ical point connects a liquid state of small strings to a 
phase with an ordered state made up of infinitely long 
strings. Thus it is natural that the transition point is 
described by a liquid state of long strings. From this pic- 
ture, one might expect the same phenomenon to occur 
in (3 + 1)D spin models. In that case, one would expect 
an entire deconfined phase between a small string phase 
and a striped phase. 



where b',c',d',e' are the renormalizcd continuum coef- 
ficients corresponding to b, c,d,e. The transition is de- 
scribed by tuning d through 0. The critical theory is 
thus a Rokhsar-Kivelson type theory described by a U(l) 
gauge theory with dynamical exponent z — 2. 

The above is the mean field critical theory. In princi- 
ple, we need to consider fluctuations to get the full critical 
theory. In particular, we need to include instanton fluctu- 
ations, replacing the (V x A) 2 term with cos(V x A). Wc 
also need to include higher order terms in E, A (which 
are particularly important in the stripe (c' < 0) phase). 
The most general terms consistent with the lattice sym- 
metry are given by 

C! = -f'E 4 - g'((E ■ ni )(E ■ n 2 )(E ■ n 3 )) 2 + ... (22) 

where ni,n.2,ri3 are unit vectors along the three lattice 
directions. 

These fluctuations have an important effect on both 
phases. The instantons destroy the string condensation 
on both sides of transition. The result is that the tran- 
sition is actually a simple Z & symmetry breaking tran- 
sition between a small string phase and a stripe phase 
(Fig. 7(b)). The higher order terms £ are also impor- 
tant - determining the ultimate form of the ordering in 
the stripe phase. 

However, at the critical point, both the instantons 
and the higher order terms £ are irrelevant for e' suffi- 
ciently small. The critical theory is therefore described 
by the simple mean field Lagrangian (21). This result 
follows from the analysis in Ref. [18]. In that paper, 
the authors analyzed the z — 2 critical point (21) in 
the context of a quantum dimer model on the honey- 
comb lattice. They found that the instanton fluctuations 
were irrelevant for e' sufficiently small. Moreover, they 
found only one relevant higher order term - a cubic term 
(E-ni)(E-ri2)(E-n 3 ). The same analysis can be applied 
in our case, but the cubic term is not allowed because of 
the symmetry E — > —E. 



VI. CONCLUSION 

In this paper, we have described a mean field technique 
for quantum string (or dimer) models. The technique can 
be used to estimate phase diagrams, to analyze the low 
energy dynamics in each of the phases, and to understand 
the critical points separating them. 

The new mean field theory developed here is more pow- 
erful than traditional mean field approaches in that it is 
applicable to both string condensed phases and the usual 
symmetry breaking phases. Thus it can be used to study 
phases and phase transitions beyond Landau's symmetry 
breaking paradigm. One particularly interesting applica- 
tion is to frustrated spin systems with emergent photon- 
like excitations. 

We have demonstrated the approach with a simple ex- 
ample: the XXZ model (1) in the limit J\ ^> J xy ^> 
\J\ — Ji\- In that limit, the low energy physics of the XXZ 
model is described by a quantum string model which can 
be studied using the mean field theory. We find that 
the model is in a paramagnetic phase for large positive 
J/g, in a plaquette ordered phase for large negative J/g, 
and that the phase transition is in the 3D XY universal- 
ity class. We have also applied the mean field approach 
to the XXZ model (1) with an additional next-nearest 
neighbor coupling ^3^((ij)) The mean field the- 

ory predicts that the model is in a paramagnetic phase 
for large positive J/g, in a stripe ordered phase for large 
negative J/g, and that the phase transition is a z = 2 
deconfined quantum critical point (21). 

Given these results, it would interesting to study the 
XXZ model (1) numerically. A quantum Monte Carlo 
study could potentially access the z = 2 deconfined quan- 
tum critical point. In addition, it could resolve the dis- 
crepancy between the resonating plaquette phase pre- 
dicted by Ref. [10] and the frozen plaquette phase that 
appears to be favored by the mean field approach (Fig. 
8). 



A natural direction for future research would be to ap- 
ply the mean field approach to a (3 + 1) dimensional spin 
model. Indeed, the mean field method may be most use- 
ful in this context. While string condensed phases with 
emergent photons are always unstable in (2 + 1) dimen- 
sions, there is no such problem in (3 + 1) dimensions. 
Thus, the mean field method can be used to find entire 
phases with string condensation and emergent photons, 
in addition to deconfined quantum critical points. 
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APPENDIX A: MINIMIZING THE ENERGY 
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FIG. 14: A plot of the minimum energy a as a function of 
J/g, obtained from a variational Monte Carlo simulation on 
a 12 x 12 lattice. The (mean field) critical point separating 
the string condensed and small string phases is given by the 
value of (J/g) where a = a c m 0.84. 



In general, the mean field phase diagram should be 
computed by minimizing the ground state energy (H) 
over all choices for and then identifying the quan- 

tum phase associated with the minimum energy ^{z}- 
Ideally, the minimization of (H) should be done in an 
unbiased fashion and include all possible {z}. However, 
to simplify our numerics, we have executed a more re- 
stricted minimization. In this section, we describe this 
minimization procedure. 

In the restricted minimization procedure, we only min- 
imize (H) over the string liquid states, zy = a. We then 
check to make sure that the minimal zy — a state is 
stable to infinitesimal perturbations a — > a + 8z%j . 

To do this, we parameterize the perturbations by a + 
Szij = ae Ei t +lAi ' where £y and Ay are real - as in 
section IV. Expanding the energy to quadratic order in 
E, A, one finds 

(H) = ^ (cij,kiEijEki + dij tk iAijA k i) (Al) 

ij,kl 

where the constants Cij,kl'dij,kl are given by the equal 
time correlation functions (8). 

To check for stability one needs to check whether the 
matrices cy,fej , dij,kl are positive definite. This can be ac- 
complished most easily by going to Fourier space, where 
cy ,ki > dij t ki are diagonal. If all of the resulting eigenval- 
ues are positive, then the state *Sf a is stable. Otherwise 
it is unstable. 

If it is stable, we assume that it is the true lowest 
energy state. If it is unstable, say in the direction a — > 
ae Ei i +lAi i 1 we conclude that the system enters a new 
(symmetry breaking) phase with an ordering given by 
ae EiS+tAii . In effect, by restricting attention to local 
instabilities, we assume that the phase transitions out of 
the liquid phases are second order or weakly first order. 

This restricted minimization procedure is less power- 
ful and less reliable then a general minimization of (H) 



over all ^{ z }- Its only advantage is that it is technically 
simpler to implement. 



APPENDIX B: CALCULATION OF MEAN FIELD 
PHASE DIAGRAM 

In the following, we describe in detail how the mean 
field phase diagram for (3) was obtained. We begin with 
the transition at J/g = 0.27. The transition at J/g = 
0.27 was obtained by calculating the optimal (minimal 
energy) a for different values of J/g. Plotting a as a 
function of J/g, we find that when J/g is reduced below 
0.27, the minimum energy a becomes larger than a c = 
2-V4 = 0.8409... (see Fig. 14). Thus, at this point, 
the variational ground state enters the string condensed 
phase. 

The string condensed (liquid) phase persists for 
-0.43 < J/g < 0.27. When J/g < -0.43, we find that 
the liquid becomes unstable (see Fig. 15). In fact, two 
modes become unstable simultaneously. The two unsta- 
ble modes are of the form 

E rs = E ±Q . sr = -j^=(n ■ (s - r)e ±iQ - r (Bl) 
V3Na A 



where ±Q are the two wave vectors shown in Fig. 12, 
and nlQ. 

Because two modes become unstable simultaneously, 
there is an ambiguity in the way the system orders. The 
system could potentially order in any (real) linear com- 
bination of the two modes: E rs = £E + Q rs + E^Q rs . 

A purely quadratic analysis cannot distinguish be- 
tween this continuum of possible orderings. However, we 
expect that higher order terms will pick out a particular 
ordering. Indeed, let us consider the energy as a function 
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FIG. 15: The eigenvalue cq corresponding to the mode (Bl), 
as a function of (J/g). The liquid becomes unstable when cq 
becomes negative. From the above plot, we find an instability 
at (J/g) c ~ -0.43. 



of the complex parameter £, (H) = H((). The symme- 
try of the lattice requires that H(Q = H(-Q = H((*) 
and H(C) = H(oj() where w is a third root of unity. The 
most general form for H satisfying these constraints is 

H(() = A|C| 2 + B|C| 4 + C(( e + (Cf) + ... (B2) 

to sixth order in £. The first two terms are not sensitive to 
the phase of £ and therefore tell us nothing about which 
linear combination is favored. However, the third term 
does pick out a phase. The phase depends on whether 
C is positive or negative. If C is positive, imaginary £ is 
favored. This corresponds to a "frozen plaquette" phase - 
where a third of the plaquettes are typically occupied by 
strings, and the interstitial bonds are likely to be empty 



(Fig. 8a). On the other hand, if C is negative, real £ 
is favored. This corresponds to a "resonating plaquette" 
phase where a third of the plaquettes resonate between 
two different configurations, while the interstitial bonds 
are typically occupied (Fig. 8b). 

While it is difficult to determine which of these two 
possibilities occur using the restricted minimization pro- 
cedure, one can make a determination by implementing 
a general minimization of (H) over all ^{ z }- We have 
implemented this on a 3 x 3 lattice and the result is that 
the frozen plaquette phase is favored. Based on this re- 
sult, we believe that C is positive and the instability at 
J/g = —0.43 results in frozen plaquette order. However, 
a more complete numerical study is necessary to even 
make a definitive mean-field prediction. 

When J/g is decreased below —0.43 the system ac- 
quires (frozen) plaquette order. However, the ordering is 
weak near the transition point, and a finite amount of or- 
dering is necessary to destroy string condensation. This 
means that the string condensation persists for a finite 
interval below —0.43 - the system enters into a phase with 
simultaneous plaquette order and string condensation. 

We have not executed systematic numerics beyond 
—0.43. However, small lattice results suggest that pla- 
quette order strengthens as J/g decreases. This sug- 
gests that when J/g becomes sufficiently large and neg- 
ative, the plaquette order becomes sufficiently strong 
that string condensation can no longer coexist and is de- 
stroyed. The system then enters a phase with the same 
plaquette order but no string condensation. 

The phase diagram for the model (1) with next nearest 
neighbor coupling ^3^((ij)) was computed using 

the same technique. We will not repeat the details here 
because of their similarity to those described above. 
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